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CN ' Abstract 

Systems that involve TV identical interacting particles under quantum confinement appear through- 

O ■ 

out many areas of physics, including chemical, condensed matter, and atomic physics. In this paper, we 
present the methods of dimensional perturbation theory, a powerful set of tools that uses symmetry to 
yield simple results for studying such many-body systems. We present a detailed discussion of the dimen- 
sional continuation of the TV-particle Schrodinger equation, the spatial dimension D — > oo equilibrium 
(D°) structure, and the normal-mode (-D -1 ) structure. We use the FG matrix method to derive general, 
analytical expressions for the many-body normal-mode vibrational frequencies, and we give specific an- 



alytical results for three confined TV-body quantum systems: the TV-electron atom, TV-electron quantum 
dot, and TV-atom inhomogeneous Bose- Einstein condensate with a repulsive hardcore potential. 

1 Introduction 

During the last two decades, techniques to confine and manipulate atoms, ions, and electrons have led to the 
creation of new TV-body quantum systems of both fundamental and technological interest. The condensation 
of atomic bose gases, the confinement of atoms using optical lattices, and the creation of quantum dots in 
semiconductors are all examples of novel many-body environments that can be manipulated using external 
potentials. These systems thus provide a unique opportunity to study many-body effects over a range of 
interaction strengths as external potentials are tuned. The quantum dot is a nanostructure in which electrons 
are confined in all three dimensions by an external potential. This "artificial atom" has shown great potential 
in biological and medical applications and may form the basis for a new generation of semiconductor lasers 
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with applications in quantum computation and quantum cryptography. The atomic vapor Bosc-Einstcin 
condensate, arising from the collapse of a collection of identical, harmonically trapped bosons into the 
lowest single-particle state of the harmonic oscillator potential, provides a source of coherent matter waves 
with possible applications in atom interferometry, atom circuits, holography, precision measurements, and 
quantum computation. The confinement of atoms in single cells of optical lattices could provide the controlled 
interaction required to create a quantum logic gate. 

These new quantum confined TV-body systems have resulted in renewed interest and new demands in 
the many-body techniques of quantum physics and chemistry, originally developed to study atoms and 
molecules - TV-body systems confined by the mutual attraction of the particles themselves. Mean-field 
treatments of these TV-body systems, such as the Hartree-Fock method in atomic physics and the Gross- 
Pitaevskii approximation for Bose-Einstein condensates, do not include many-body effects, and, therefore, 
become inaccurate for systems under tight confinement or strong interaction. These new systems, which can 
have a few hundred to millions of particles, present serious challenges for existing many-body methods, most 
of which were developed with small systems in mind. 

In this paper we offer the methods of many-body dimensional perturbation theory for the study of 
large TV-body systems under quantum confinement. We apply dimensional perturbation theory (DPT) to 
many-body quantum confined systems from chemical, condensed matter, and atomic physics. The first 
application of many-body DPT (Section 5), and the impetus for this paper, originated in Ref.[l], in which 
Loeser introduced low-order many-body dimensional perturbation methods to the TV-electron atom. In this 
instrumental paper, Loeser obtains low-order, analytical expressions for the ground-state energy of neutral 
atoms. For Z = 1 to 127, the numerical results compare well to Hartree-Fock energies with a correlation 
correction. In Sections 2-4 of this paper, we discuss the formalism and general theory behind Loeser's results 
for the TV-electron atom, in which the quantum confinement of the TV electrons is supplied by the Coulomb 
attraction of a nucleus. The second application is the quantum dot (Section 6), an atom-like many-body 
system from condensed matter physics, where the confinement of the TV electrons is supplied by an external, 
isotropic trapping potential. The final system (Section 7) is TV identical hard spheres in an isotropic trap. 
Although our results are more general, this hard-sphere model is appropriate for describing a Bose-Einstein 
condensate of a trapped alkali- metal gas, a system that, like the quantum dot, bridges the areas of atomic 
and condensed matter physics. 

In many-body dimensional perturbation theory, the coordinate vectors of the TV particles are allowed to 
have D Cartesian components. The singular limit D — > oo is taken as the unperturbed system, and l/D be- 
comes the perturbation parameter. The Schrodinger equation is written in terms of a similarity transformed 
wave function that removes first-order derivatives from the Laplacian and introduces a centrifugal-like poten- 
tial containing all of the explicit dimension dependence of the Laplacian. The explicit dimension-dependence 
of the centrifugal-like potential is quadratic, and in order to regularize the large-dimension limit of the Hamil- 
tonian (the D — > oo limit or the leading term in the l/D expansion) we choose a scaling of the length and 
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energy that is also quadratic in D. Then as D — > oo the second derivative parts of the kinetic energy 
vanish and the particles become localized in the bottom of an effective potential defined by a centrifugal-like 
contribution from the kinetic energy and contributions from the other potential energies (i.e., the confine- 
ment and interaction potentials). In the case of attractive interparticle forces, the repulsive centrifugal- like 
term stabilizes the large- D configuration against collapse. The first-order quantum correction corresponds 
to normal-mode vibrations about the large- D effective potential minimum. We find the normal-mode fre- 
quencies using the FG matrix technique of molecular physics [2]. Higher orders can be calculated using a 
matrix method developed specifically for dimensional perturbation theory [4]. 

This paper is organized as follows. In Section 2 we give the unsealed Schrodinger equation, and in Section 
3 we discuss scaling methods and the zeroth-order (D — > oo) term of the 1/D expansion. In Section 4 we 
discuss the 1/D expansion and utilize the FG matrix method, and we derive analytical expressions for the 
normal-mode frequencies of vibration about the infinite- Z) symmetric minimum, from which we derive the 
first-order energy approximation. In Sees. 5-7, we apply these methods to the TV-electron atom, TV-electron 
quantum dot, and iV-atom Bose-Einstein condensate. 



2 D-dimensional TV-body Schrodinger equation 



For an TV-body system of particles confined by a spherically symmetric potential and interacting via a 
common two-body potential gij , the Schrodinger equation in D-dimensional Cartesian coordinates is 



N N-l N 

i—1 i—1 j— i-\-l 



where 



D 



hi 



V c< 




and = V int ^y^„ = i ^ Xiv ~ x ^f^j 



(1) 



(2) 



(3) 



arc the single-particle Hamiltonian and the two-body interaction potential, respectively. The operator H 
is the D-dimensional Hamiltonian, and x iv is the v th Cartesian component of the i th particle. V conf is the 
confining potential. For the iV-electron atom, V conf is the Coulomb attraction between the electrons and the 
nucleus, while for the TV-electron quantum dot and TV-atom hard-sphere problem we model the confinement 
as a harmonic trapping potential. The two-body interaction potential V± nt is Coulombic in the first two 
systems and a hard sphere in the third. 



3 



2.1 Transformation of the Laplacian 

Restricting our attention to spherically symmetric (L = 0) states of the many-body system, we transform 
from Cartesian to internal coordinates. A convenient internal coordinate system for confined systems is 



n 



\i x i" ( l <i< N ) and lij = cos(%) = XivXj^J /ntj (1 < i < j < N), (4) 

which are the D-dimensional scalar radii j-j of the N particles from the center of the confining potential and 
the cosines 7^ of the N(N — l)/2 angles between the radial vectors. 

Now for a function \I> dependent on only two functions r(x) and -f(x) one can write 

d 2 t>(r(x),j(x)) _ d 2 r d<f d 2 1 d^> fdr\ 2 d 2 t> ( dj\ 2 d 2 ^ dj dr d 2 ^ 



dx 2 dx 2 dr dx 2 d^y \dx J dr 2 \dx J dj 2 dx dx drdj 

Generalizing this, when operating on the state $?(ri(xi v ), Jn(xi v ) . . . 7ik{xi v ) . . . JiN{xiv)) (where k ^ i and 
v = !,...,£))), Vf can be written in terms of the internal coordinates of Eq. (4) as 



/ dri \ ( djij \ d 2 



\dx iv ) \dx iv ) drtdfij^' 
The relevant derivatives of the internal coordinates are 



2-— 1_ T2" q 2 = I2 I 3 ~ ■- 2 — — -7*j J > (8) 



ftc* V rf J dxiS rf \ rf n r 3 

which lead to the effective S'-wave Laplacian in internal coordinates: 



I) 1 d , \ -> D—l x d , 



2.2 Removal of first-order derivatives 

Next we wish to find a transformation of the Hamiltonian that removes the first-order derivatives from the 
Laplacian (Eq. (9)) so that the kinetic energy operator is reduced to a sum of terms of two kinds, namely, 
a second-order derivative term and a repulsive centrifugal-like term, which when attractive interparticle 
potentials are present, stabilizes the system against collapse in the large- D limit. When this is done the 
zeroth and first orders of the dimensional (1/D) expansion of the Hamiltonian become exactly soluble for 
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any value of N. In the D — > oo limit, the second derivative terms drop out, resulting in a static problem 
at zeroth order, while first order corrections correspond to simple harmonic normal-mode oscillations about 
the infinite-dimensional structure. 

In Rcf. [3] , Avery et al. considered the problem of performing a similarity transformation of the wave 
function ^ and operators O: 

$ = X~ 1 *, and d = X - 1 d X , (10) 
where the transforming function, with adjustable parameters a and (3, is of the form: 

X=(nr 2 ...r N )- a T-^ 2 . (11) 

Here T is the Gramian determinant, the determinant of the matrix whose elements are jij (see Appendix 
D). One of the cases considered by Avery et ai.[3] for a and (3 (ai = (D — l)/2 and (3i = (D — N — l)/2) 
causes the weight function for matrix elements, W = J\ 2 , to equal unity, where J is the Jacobian of the 
transformation to internal coordinates: 

J=(r 1 r 2 ...r N ) D - 1 ^ D - N - 1 ^ 2 . (12) 

As we discuss in another publication [5], this scenario has the advantage of making matrix elements easier to 
calculate using DPT and makes the physical interpretation of the large-dimension normal-mode structure 
more transparent. However, the results in this paper are easier to derive if we use a = (3 = (D — l)/2, which 
removes the first-order derivatives from the Laplacian while giving the same results as ol\ and (3\ through 
first order in 1/D. Carrying out the transformation of the Hamiltonian and wave function of Eq. (1) via 
Eq. (10) with a = (3 = (D — l)/2 in Eq. (11), the Schrodingcr equation becomes[l]: 

(T + V)<I> = E<£ (13) 

^ 92 y^ ft 2 (7jfc-7ij7ifc) d 2 h 2 (D — 1){D — 2N — 1) £W \ 
2m, dr 2 2 mi r 2 d^d^k 8m<r? T [ ' 

N N-l N 

v = J2 VoMn) + E E v^inj). (15) 

The Gramian matrix whose determinant is r« is the i th principal minor formed by deleting from T the row 
and column corresponding to the i th particle. The quantity = ^Jrf + r 2 — 2rirjjij is the interparticle 
separation. The similarity transformed Hamiltonian for the energy eigenstate $ is x _1 i?x, where H = 
(T + V). 

We remark that for the hard-sphere system in Sec. 7, we find it expedient to choose a dimensional 
continuation of Vi nt that contains explicit dimension dependence, which is not expressed in the equation 
above. However, the general discussion to follow in Sees. 3-4 holds for the hardsphere system with only 
slight modification to the first-order energy approximation. This is discussed further in Sec. 7. 
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3 Infinite-/} analysis: Leading order energy 

To begin the perturbation analysis we regularize the large-dimension limit of the Schrodingcr equation by 
defining dimcnsionally scaled variables: 

h = n/n{D) , E = k(D)E and H = k{D)H (16) 

with dimension-dependent scale factor k(D). From Eq. (14) one can see that the kinetic energy T scales in 
the same way as 1/r 2 , so the scaled version of Eq. (13) becomes 

H<P = ( — 5— f + V ) $ = E$, (17) 



AD) 

where barred quantities simply indicate that the variables are now in scaled units. Because of the quadratic 
D dependence in the centrifugal-like term in T of Eq. (14), we conclude that the scale factor k(D) must also 
be quadratic in D, otherwise the D — > oo limit of the Hamiltonian would not be finite. The factor of k(D) 
in the denominator of Eq. (17) acts as an effective mass that increases with D, causing the derivative terms 
to become suppressed while leaving behind a centrifugal-like term in an effective potential, 



N , 2 p(i) x JV-1 AT 

= E ( ^2 - jr + Wf*)J + E E ( 18 ) 



in which the particles become frozen at large D. In the D — » oo limit, the excited states have collapsed onto 
the ground state, which is found at the minimum of V el ± . 

We assume a totally symmetric minimum characterized by the equality of all radii and angle cosines of 
the particles when D — > oo, i.e., 

f i =f 00 (1 < * < TV), = Too (1 < * < j < N). (19) 

Since each particle radius and angle cosine is equivalent, we can take derivatives with respect to an arbitrary 
fi and 7y in the minimization procedure. Then evaluating all f, and 7^ at the infinite-D radius and angle 
cosine, and 700, respectively, we find that and 700 satisfy 



dV ei 



dn 
dV ell 



djij 



= (20) 
= 0, (21) 



where the 00 subscript means to evaluate all rj at and all jij at 7^. In scaled units the zeroth-order 
(D — » 00) approximation for the energy becomes 



E ao = V e 



eff 



= Wfoo,7oo). ( 22 ) 



In this leading order approximation, the centrifugal- like term that appears in V S i± , even for the ground state, 
is a zero-point energy contribution required by the minimum uncertainty principle [6]. 
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4 Normal- mode analysis: 1/D first-order quantum energy correc- 
tion 



At zeroth-order, the particles can be viewed as frozen in a completely symmetric, high-D configuration 
or simplex, which is somewhat analogous to the Lewis structure in atomic physics terminology. Likewise, 
the first-order 1/D correction can be viewed as small oscillations of this structure, analogous to Langmuir 
oscillations. Solving Eqs. (20) and (21) for foo and gives the infinite- -D structure and zeroth-order energy 
and provides the starting point for the 1/D expansion. To obtain the 1/D quantum correction to the energy 
for large but finite values of D, we expand about the minimum of the D — > 00 effective potential. We first 
define a position vector, consisting of all N(N + l)/2 internal coordinates, whose transpose is: 



V = {n,r 2 , ■ ■ • ,r N , 712, 713, ■ • - ,7jv-i,jv)- 
We then make the following substitutions for all radii and angle cosines: 

f i = f 00 + jvy. 

Hi =7=o +£ 1/2 7i,-, 



(23) 

(24) 
(25) 



where S = 1/D is the expansion parameter, and we define a displacement vector consisting of the internal 
displacement coordinates [primed in Eqs. (24) and (25)]: 



r.'T 



y ={r' 1 ,r' 2 ,...y N ,i 12 ,i 13 ,...,i N _ hN ). (26) 

We may then obtain a power series in 5 1 / 2 of the effective potential about the D — > 00 symmetric minimum: 

d 2 V eii 



V elf = V. H 


+ < 5l/2 E 


\dV eM 

[ Qy'a _ 






OO ( 

(7=1 




00 



dy' a dy' v 



O 



(^ 3/2 ) , 



(27) 



where 



P = iV(7V+ l)/2 (28) 

is the number of internal coordinates. The first term in the power series (Eq. (27)) is simply the zeroth- 
order energy (Eq. (22)). The second term is zero since we are expanding about the minimum of the effective 
potential; the system is said to be in equilibrium since the forces acting on the system vanish [Eqs. (20) 
and (21)]. The third term defines the elements of the Hessian matrix[7] F of Eq. (30) below. The derivative 
terms in the kinetic energy are taken into account by a similar series expansion, beginning with a first-order 
term that is bilinear in d/dy', i.e., 

p p 

T = 2 5 E E Gau^dy, + O (> 2 ) , (29) 

(7=1 V=\ 

where T is the derivative portion of the kinetic energy T [see Eq. (14)]. Thus, obtaining the first-order 
energy correction is reduced to a harmonic problem, which is solved by obtaining the normal modes of the 
system. 
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We use the Wilson FG matrix method [2] to obtain the normal- mode vibrations and, thereby, the first- 
order energy correction. It follows from Eqs. (27) and (29) that G and F, both constant matrices, are 
defined in the first-order S = 1/D Hamiltonian as follows: 



--d v ,+ Gdy, + -y' T Fy' 



(30) 



After the Schrodinger equation (13) has been dimensionally scaled, the second-order derivative terms are of 
order 6, and by comparing these terms with the first part of Hi, the elements of the kinetic-energy matrix 
G are easily determined. The elements of the Hessian matrix[7], F av = [d 2 V e ff/dy' a y' u ]oo, on the other hand, 
require a bit more effort to obtain, as we will see in detail in later sections. 

We include a derivation of the FG matrix method in Appendix A, but we state here the main result of 
the method, which consists of finding the roots of the following characteristic polynomial in A: 



det(AI - GF) = 0. 



(31) 



Depending on the number of particles, the number of roots A - there are P = N(N + l)/2 roots - is 
potentially quite large. However, as we will see shortly, there is a high degree of degeneracy due to the total 
symmetry of the infinite- D Lewis structure. In fact, there are only five distinct roots, A M , where [i runs over 
, + , 1 , 1 + , and 2. And as we conclude in Eq. (60), the energy through first-order can be written in 
terms of the distinct normal-mode vibrational frequencies, which are related to the roots A M of GF by 



\ _ -2 



(32) 



as can be seen from Eq. (A5) in Appendix A. In the next subsections, we show explicitly how to find 
analytical expressions for the roots of GF. 



4.1 Indical structure of F, G, and GF matrices 

The F, G, and GF matrices, which we generically denote by Q, are P x P matrices with the same indical 
structure as yy T : 
( 



yy 

















r\r 2 ■ 


nrjv 


ri7i2 


ri7i3 


ri7Ar_iAr 




T2T2 ■ 




^2712 


^2713 


T21N-1N 


f N fl 






rjv7i2 




rNlN-lN 


712^1 


712^2 ' 


7i2^v 


712712 


712713 


l\2lN~lN 


713^1 


713^2 ' 


712^1 


713712 


713713 


7137N-1N 


jN-lNn 




• In-inTn 


7W-17V712 


7N-1AT713 • 


• IN-IN^N-IN 



(33) 



where y is defined by Eq. (23). The indical structure of this matrix suggests a convenient shorthand for 
referencing the elements of the Q matrices. The upper left block of Eq. (33) is an (N x N) matrix with 
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elements associated with (ft, fj); hence we use the subscript (i, j) to refer to these elements. The upper right 
block is an (N x N(N — l)/2) matrix with elements associated with (^,7^); hence, we use the subscript 
(i,jk) to refer to these elements. The lower left block is an (N(N—l)/2 x N) matrix with elements associated 
with (-fij,fk); hence, we use the subscript (ij, k) to refer to these elements. Finally, the lower right block is 
an (N(N — l)/2 x N(N — l)/2) matrix with elements associated with (7y,7fcj); hence, we use the subscript 
(ij, kl) to refer to these elements. 

4.2 Symmetry of the Q matrices 

As the number of particles N increases, diagonalizing the P x P GF matrix (where P = N(N + l)/2) 
becomes, prima facie, a daunting task. However, one of the advantages of dimensional perturbation theory is 
the simplifications that occur in the large-dimension limit. In particular, since we are dealing with identical 
particles in a totally symmetric configuration (the Lewis structure) in which all the particles are equivalent, 
the Q matrices display a high degree of symmetry with many identical elements. Specifically, 

Qa 

Q b (i ? j) and (i' ± f) 

Q c (i ± j) and (i' ± f) 

Q d (i j ? k) and (*' ? f ? k') 

Q e (i j) and (i 1 ? f) (34) 

Q f (i * j + k) and (*' ? f ? k') 

Q g (i j) and (i' f) 

Q h (i ? j k) and (»' ? f k') 

Q t (i + 3 + k + V) and (i' ^ f ^ k' ^ I'). 

Note the indices in the relationships above run over all particles (1,2, ... ,N) with the exceptions noted in 
the far right column. For example, Qij — Qi\j> = Qb, where (i ^ j) and (i' ^ j'), means that all off-diagonal 
elements of the upper left block (the pure radial block) of Q are equal to the same constant Qb. Similarly, 
Qij,ki — Qi'j',k'l' = Ql, where (i 7^ j 7^ k ^ I) and (i' 7^ j' 7^ k' 7^ I'), means that any elements of Q in 
the lower right block (the pure angular block) that do not have a repeated index are all equal to the same 
constant Q L . We should remark here that G and F are also symmetric matrices (G T = G and F T = F); 
however, while GF does display the high degree of symmetry of Eq. (34), it is not a symmetric matrix. 

4.3 Q matrices in terms of simple submatrices 

The symmetry of the Q matrices (F, G, and GF) described in Eq. (34) allows us to write these matrices in 
terms of six simple submatrices. We first define the number of 7^ coordinates to be 

M = N(N -l)/2, (35) 





— Qi' ,i' 




Qi.j 


Qi f ,j f 




Qij,i 


— Qi'j'^i' 




Qjk,i 


Qj'k',i' 




Qi.ij 


Qi'.i'j' 




Qi.jk 


= Qi'j'k' 




Qij,ij 


— Qi'j'.i'j' 




Qij,jk 


— Qi'j'J'k' 




Qij,kl 


— Qi'j'.k'U 
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and let I at be an N x A" identity matrix, Im an M x M identity matrix, Jn an N x N matrix of ones and 
3m an M x M matrix of ones. Further, we let R be an N x M matrix 1 such that Rijk = &ij + <5ifc, Jjvm 
be an JV x M matrix of ones, and 3 NM = Jmn- 
We then write the Q matrices as 

Q-f Ql *). <»> 

V Q 3 Q 4 J 

where the block Qi has dimension (N x N), block Q 2 has dimension (N x M), block Q 3 has dimension 
(M x iV), and block Q4 has dimension (M x M). Now, as we show in Appendix B, Eq. (34) allows us to 
write the following: 

Qi = (Q a - Q h )l N + Q h J N (37) 

Q 2 = {Qe-Qf)~R + QfJNM (38) 
Q 3 = (Qc - Qd)R T + Qd^NM (39) 

Q 4 = (Q g -2Q h + Q l ,)i M + (Qh-Q,)n T li + QJM. (40) 

In particular, letting Q = GF, the matrix that must be diagonalized, Eq. (36) becomes 

GF= ( (a-b)I N + bJ N (e-f)R + fJ N M jjr 

y (c-d)R T + dJ M N (g -2h + l)I m + {h- l)TL t H + ii M 

where we have used the following abbreviations: 



a 




(GF) a 


— G a F a 


b 




(GF) b 


= G a F b 


c 




(GF) C 


= G g F e + {N- 2)G h F e + {N- 2)G g F f 


d 




(GF) d 


= G g F f + 2G h F e + 2(N - 3)G h F f 


e 




(GF)e 


= G a F e 


f 




(GF) f 


= G a F f 
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(GF)g 


= G g F g + 2(N-2)G h F f 


h 




(GF) h 


= G g F h + G h F e + (N- 2)G h F h + (N - 


1 




(GF) L 


= G g F b + 4G h F h + 2(N - A)G h F t . 



(42) 



The right-hand sides of Eq. (42), the GF matrices expressed in terms of the F and G matrix elements, are 
derived in Appendix B. 

1 In graph theoretic terminology R is called a vertex-edge matrix (see Appendix B). 
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4.4 Normal mode frequencies and first-order energy 

Having obtained the GF matrix of Eqs. (41) and (42), we may then, according to Eq. (31), find the 
eigenvalues A of GF by solving 

det(E) = 0, (43) 

where we have defined E as 

/ AIjv \ 

E= -GF. (44) 

\ Mm J 

To find an analytical expression for det(E) we multiply E by three matrices: X, Y and Z, such that 

/ ±1 N + ^J N \ 

XYEZ = " v , (45) 

y • • • vi M J 

where det(X) = det(Y) = det(Z) = 1 so that det(E) = det(XYEZ). In Appendix C we construct the 
matrices X, Y and Z needed to transform E of Eq. (44) to XYEZ of Eq. (45). We find 

v = A - g + 2h - i (46) 

* = {X-a + b)v + {N-2)(t-h)(X-a + b) + (N-2)(d-c)(e-f) (47) 
u = -kv- (h- a + b) + (d- c)(e- /) 

+ [f(2(d - c) - Nd) -i{\-a-{N - 1)6) + 4b(h - b) - 2d(e - /)] . (48) 

The determinant of XYEZ, which does not depend on the submatrix in the lower left block, is found 
from Appendix B [see Eq. (B30)] to be 

det(E) = dct(XYEZ) = ^^(t + Nu)v M ~ N . (49) 

There are five distinct roots, which naturally arise from the terms in Eq. (49) with t, u, and v given by Eqs. 
(46), (47) and (48). From Eq. (49), det(E) factors into three terms, resulting in three equations for the five 
distinct roots: 

t N - x = (50) 
(t + Nu) = (51) 

y M-N = v N(N-3)/2 = Q; (52) 

From Eqs. (46) and (47), one can see that t is quadratic in A. Hence, from Eq. (50) there are two roots, 
which we denote by 1~ and 1+, with multiplicity d 1 - = c?]+ = N — 1 given by 

t = 0. (53) 

Similarly, the (t + Nu) term in Eq. (49) is quadratic in A and, hence, there are two roots, which we denote 
by 0~ and + , with multiplicity d - = d + = 1. Finally, v is linear in A; hence, there is one root, which we 
denote by 2, with multiplicity d 2 = N(N — 3)/2 from Eq. (52) given by 

v = 0. (54) 
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From Eqs. (54) and (46), the root of multiplicity N(N - 3)/2 is 

\ 2 = g — 2h + l. 



(55) 



From Eqs. (53) and (47), the two roots with multiplicity (N — 1) are 

Ai± = r/i ± v^i 2 - Ai, 

where 



(56) 



(57) 



m = \ [a - b + g + {N - A)h - {N - 3)t] 

Ai = (N - 2)(c - d)(e - /) + (a - b) [g + (N - 4)h - (N - 3)i] . 
From the i + TVu term in Eq. (49) along with Eqs. (47) and (48), the two roots with unit multiplicity are 



A ± 



A fi 



where 



Vo = 



a — (N — l)b + g + 2(N - 2)h + (JV ~ 2) 2 (JV ~ 3) t 



A = (a - (JV - 1)6) \g + 2(N - 2)h - ( N ~ 2 W-V L ] _ JY_2 (2c + (JV _ 2)rf)(2g + (JV _ 2)/) 



(58) 



(59) 



The five distinct roots A^, where /x runs over , + , 1 , 1 + , and 2, are equal to the square of the 
normal- mode vibrational frequencies (i.e., Eq. (32)), and the energy through first-order in 6 = 1/D is then 



E 



E^ + SE + 0(8 2 ) 



OO ^ 



(60) 



^={0 ± . 1±, >v=o 
2 } 



where the n M are the vibrational quantum numbers of the normal modes of the same frequency u)^ (as such, 
n M counts the number of nodes in a given normal mode). The quantity d^.„^ is the occupancy of the manifold 
of normal modes with vibrational quantum number n M and normal mode frequency u> M , i.e. it is the number 
of normal modes with the same frequency lj^ and the same number of quanta n M . The total occupancy of 
the normal modes with frequency is equal to the multiplicity of the root A p (see the discussion after Eq. 
(52)), i.e. 

OO 

n M =0 

where is the multiplicity of the /i th root. Note that although the equation in Ref. [1] for the energy 
through Langmuir order is the same as Eq. (60), it is expressed a little differently (See Ref. [8]). 



5 iV-electron atom 

The ./V-electron atom was previously discussed using low-order many-body dimensional perturbation theory 
by Loescr[l, 20]. In this section, we give details necessary to derive the results in Ref. [1], to which we direct 
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the interested reader for specific numerical results. The confining potential is provided by the Coulomb 
attraction of the electrons to the nucleus: 



Vconf {n) 



z 

1 

r,. 



and the interaction potential is the Coulomb repulsion of the electrons: 



f + r 2 



(62) 



(63) 



Substitution of the following charge/dimensionally scaled variables into the similarity transformed Schrodingcr 
equation (13) in atomic units (K = m = 1): 



Z Ji, E = ^E, and Q = (D — l)(D - 2N - l)/4, 



n .„ „ z2 ^ 7 ^ v _ ., v _ „. (64) 

places all of the dimension dependence in the second-derivative parts of the kinetic energy and gives the 
following scaled equation: 

(Y« + U {k) + y (A) ) $ = E$, (65) 



where 



N 



2ftE 



N 



d 2 

^ + EE 



7jk 



d 2 



-IE 



i r« 



(66) 



N AT-1 AT 

= -E-+zE E 



i=i ■ 1 - i=i j=*+i \Jr 2 +r 2 - 2r i r j7ij 

The superscript A refers to the iV-electron atom system (to distinguish these quantities from those in the 
systems to follow). From this equation, it can be seen that if D, and hence fi, become infinitely large, the 
differential part of the kinetic energy (U^ being the centrifugal-like part of the kinetic energy) will 
drop out of the similarity transformed Hamiltonian. In effect it is as though the particles become infinitely 
heavy and, because of this, any function from the basis set of all delta functions in configuration space is 
an eigenfunction of the Hamiltonian. The energy of such an eigenfunction is just the value of the effective 
potential at that specific point in configuration space which that delta function selects, and, as mentioned 
in Section 3, this point is the minimum of the infinite- 1? effective potential: 



m } = v. 



eff 



(67) 



where the effective potential v}*} corresponds to Eq. (18) in scaled units. 



5.1 iV-electron atom infinite-/} analysis 

We have chosen a particularly simple infinite- D configuration: the totally symmetric configuration, which is 
characterized by the equality of all fj and jij. The effective potential is an extremum when its derivatives 
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with respect to all r, and 7^ are zero at and 700. Using the conditions imposed in Eqs. (20) and (21), 
we find: 



1- N ^ 



1 l + (jV-2) 7co 1 

"r*, (1 + (N — l) 7oo )(l - 700) 23/2ZVT^; 
1 (2 + (TV — 2)700)700 1 1 



= 0, 

0, (68) 



T* (l + (^-l)7co) 2 (l-7oo) 2 foo 23/2Z(l - 7oo) 3 / 2 

where we have used the infinite- £) symmetric-minimum Gramian results given in Eqs. (D12) of Appendix 
D. 

One can eliminate foo from the above equations and then 700 is the negative solution of smallest magnitude 
of the following equation: 

8Z 2 7 ^ [2 + (N — 2) 7co ] 2 + 7co - 1 = 0. (69) 

For Z = N it has been verified that this extremum is a minimum. If we consider Z to be a parameter 
and we decrease this parameter, the system becomes more and more unstable and at some point the above 
extremum is no longer a minimum. This critical value Zi oca i can be obtained via Eq. (69) from the negative 
real solution of the following equation: 

(4 - 6 7c o)(l - 7=o) 3 + iV7oo(22 - 9 7c o)(l - 7oo) 2 + 207V 2 7^(1 - 7oo) + 57V 3 7 ^ = 0, (70) 

which is simply the condition under which the smallest eigenvalue of the Hcssian[7] changes sign. Zi oca i is 
always less than N for neutral atoms, and so the totally symmetric state is indeed in an energetic (local) 
minimum for neutral atoms. 

The infinite- D radius and energy can be written in terms of the solution 700 of the quartic equation (69) 
as follows: 

f« = [l + (iV-l)7oo]- 2 (71) 

E W 1 

= + {N - l) 7 co] 3 [l + (N- 2) 7 oo]/(l - 7=0). (72) 

We should remark here on the fact that we are considering a minimum in the effective potential, which allows 
us to eliminate Z from the above expressions. Specifically, we have employed the conditions (20) and (21) 
at the minimum that 

svW f)TT( k ) 

(73) 



dji 3 00 <97j. 
which leads to 

1 _ 2 3 / 2 7 oo[2 + (7V-2)7 



VT^Z <74 » 
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5.2 iV-electron atom normal modes 

We can determine the elements of G by comparing the differential terms in Eq. (30) with TW in Eq. (66) 
expanded to first order in 1/D as in Eq. (29). The non-zero elements of the G matrix are: 

G a = 1 

G 9 = 2ig^ = 2[l + (iV-l) 7co ] 4 (l+ 7co )(l- 7co ) (75) 
G h = 7 ^~ 7 r } = [l + (A-l) 7co f 7co (l- 7co ), 

where the matrix elements have been evaluated at the infinite-!? symmetric minimum and we have used the 
notation in Eq. (34). The factor of 2 in G g comes from the fact that each (ij, ij) term appears twice in 
of Eq. (66) (e.g., (21, 21) as well as (12, 12)), but it is only counted once in the G matrix (e.g., only (12, 12) 
is counted). In a similar manner to the G-matrix elements, the non-zero F-matrix elements are: 

[1 + (N - l) 7 co] 6 (2 + 2(5N - 7) 7 oo + (5JV - 6)(N - 3) 7 £, - 3{N - 1)(N - 2) 7 3 



F a 
F b 
F, 



(1 - 7co) 2 V 2 
[1 + (N — l) 7co ] 6 /- 7 oo(3 - 7co)[2 + (N - 2) 7oo ] 



(1 - 7oo) 2 V 2 
[l + (N- l) 7co ] 4 /- 7co [2 + 3(A - 2) 7co ] 



2 



(1 - 7oo) 

[1 + (N - 1) 7 , 
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Ff = 1 n / 2 J (»£) (76) 



F, 



(l-7co) 2 

[1 + (N — l) 7co ] /4 + 6(2N - 5) 7co + (12iV 2 - 53iV + 64) 7 ^ + (4iV 3 - 23iV 2 + 49N - 38) 7 ^ 



s ~ n _ „, \3 



(l-7oo) 3 V 2 
Fh = [1 + (N l)7oo] 7co (3 + {5N u)7oo + {m2 _ m + n)7 ^ ) 



(l-7oo) 3 

F = [1 tr:^°° J ^l.2 + .-V-2)-,). 

(7 7 ooJ 



where we have used Eqs. (D13) from Appendix D to evaluate the Gramian second-order derivatives. Again 
we have used the fact that we are considering a minimum in the effective potential, which allows us to use 
Eq. (74) to eliminate Z from the above expressions. Pulling out the common term among the GF matrix 
elements, we define the matrix (GF)' as: 

(GF)' = rGF, (77) 



where 

2(1 -7co) 2 



f fi. (78) 



- [1 + (N - l) 7 oc] 6 

Using the above equations along with Eqs. (42) we find the following for the elements of (GF)': 



a = 2 + 2(57V-7) 7oo + (57V-6)(AT-3)7i-3(Ar-l)(7V-2)7^ 

b = -7co(3-7co)[2 + (A-2) 7co ] 

c = -7co(l-7co)[2 + (3A-4) 7co ][2+(iV-2) 7co ][l + (iV-l) 7co ] 2 

d = 2 7 2 (l- 7co )[2 + (A-2) 7co ][l + (A-l) 7co ] 2 
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e = - 7 co[2 + 3(iV-2)7 co ][l + (iV-l) 7co ]- 2 (79) 
/ = 4 7 i[l + (iV-l)7co]- 2 

g = 2[4 + 2(47V - ll) 7oo + (47V 2 - UN + 24) 7 ^ + 3(N - 2) 7 ^] 
/i = -7oo[8 + 6(iV-3)7oo-(3iV-10hi] 
t = 87^(1-700). 

The eigenvalues of (GF)' will be equal to the eigenvalues of GF multiplied by the common factor r in Eq. 
(78). 

Using the matrix elements of Eq. (79), we obtain the two multiplicity- 1 modes from Eqs. (58) and (59): 



A$ = 5(1 - 7oo) 2 + 2iV 7 co(l - 7oo) + N 2 jL ± 



3(1 - 7co) 4 (3 + 8 7co ) - 12iV 7co (l - 7co ) 3 (l - 3 7co )~ 



\ 2N^l(l - 7co ) 2 (l - 6700) + 47V 3 7 i(l - 700) + NS. 



(80) 



For the modes with multiplicity N — 1 (A^±) and multiplicity N(N — 3)/2 (X^ 1 ) we find from Eqs. (55), 
(56), and (57): 



(Ah 



A$ = 5(l- 7 co) 2 + 9iV7oo(l-7co) + ^ 2 7L± 4 



3(1 - 7co) 4 (3 + 8 7co ) - 6iV 7co (l - 7co ) 3 (l - 8 7co )- 
\J 2iV 2 7i(l-7co) 2 (4-97co) + 3iV 3 7i(l-7co) + !iV 4 7 4 ' 

(81) 
(82) 



A^ A) = 4(1 - 7co )(2 - 5 7co ) + 27V 7co (8 - ll 7oo ) + 87V 2 7 ^. 
With r given in Eq. (78), from Eq. (32) the normal-mode frequencies u> are related to these A's by 



A(a) 



(83) 



Having obtained the normal-mode frequencies from Eqs. (80), (81), (82), and (83), the energy through first 
order is given by Eqs. (60) and (72). 



6 N- electron Quantum Dot 

We follow the ./V-electron atom with an analogous many-electron system, the quantum dot, or, as it is 
sometimes called, the artificial atom[12]. Quantum dots are nanostructures in which a controllable number of 
electrons are attracted to a central location. But instead of the Coulomb attraction of a nucleus, the quantum 
dot electrons are attracted to the center of an external trapping potential. Another difference between 
quantum dots and atoms is their size: quantum dots are typically much larger than atoms. Furthermore, 
the quantum dot electrons typically interact in some medium such as a semiconductor; thus, we use the 
effective-mass approximation where the electrons, each with mass to*, move in a medium with dielectric 
constant e. McKinney and Watson have used dimensional perturbation theory to solve for the two-electron 
quantum dot spectrum[13]. In their paper, they outline how one could use the many-body techniques 
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previously introduced by Loeser[l, 20] and detailed in this paper, to apply dimensional perturbation theory 
to a quantum dot with an arbitrary number of electrons. We now provide the details sketched in Rcf. [13] . 

The electrons, each of effective mass to*, are confined by a spherical harmonic trap with trapping fre- 
quency u> ho : 

V cont (n) = l -m*ul n\ (84) 

and we take the interelectron potential to be 

v int ( rij ) = - e2 == , (85) 

where e is the electric charge and e is the dielectric constant. Unlike the TV-electron atom problem, where we 
used dimensionally scaled atomic units, for the quantum dot we use dimensionally scaled harmonic oscillator 
units. We can transfer all of the explicit dimension dependence to the differ ential part of the kinetic energy 
and regularize the large-dimension limit by substituting the following dimensionally scaled variables into the 
similarity transformed Schrodingcr equation (13): 

E = -^—E, and Q = (£>-!)(£> -2N-l)/4, (86) 



where the dimensionally scaled harmonic oscillator length and trap frequency are, respectively, 



That is, we obtain 



where 



ko = \j and Col = ^ho- (87) 

TO Who 



(T (QD) +[/ (QD) +y (QD) )$ = j g $) (88) 



1 N 
90 ^ 



2VL 

i=l 



+ EE 



i N i rW 



i=i 

N , , N-l N 



^ D) -E^1 + 7E E 



and where 



2 ^i=i j=i+i Jrf + r 2 - Znrjjij 



eh 2 



6=7— and a* = (90) 

are the coupling constant and the effective bohr radius, respectively. From Eq. (89), it can be seen that if 
D, and hence f2, becomes infinitely large, the differential part of the kinetic energy will drop out of the 
Hamiltonian, just as for the TV-electron atom. The particles behave as though they become infinitely heavy 
and the infinite-D energy becomes the value of the effective potential, 

v m = u m + v m t (9 i) 
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at its minimum: 



£(« D ) = V 



(QD) 
eff 



(92) 



6.1 iV-electron Quantum Dot infinite-/} analysis 



As before, we choose the totally symmetric configuration for which all fj and 7^ are equal to some and 
7oo. The effective potential is an extremum when its derivatives with respect to all fj and 7^- are zero at 
this foo and 700. Using Eq. (91) in Eqs. (20) and (21), we find: 



1 



1 + (N - 2) 7o 



T* (l + (7V-l)7oc)(l-7oc) 

1 (2+(A-2) 7co ) 7co 



+ 



TV- 1 



00 

1 1 



= 0, 



0. 



(93) 
(94) 



f^(l + (iV-l)7oo) 2 (l-7oo) 2 23/^(1 -7co) 3 / 2 

where we have used the infinite- 1? symmetric minimum Gramian results given in Eqs. (D12) of Appendix 
D. 

One can eliminate from the above equations and then, just as for the iV-electron atom, 7^ is the 
negative solution of smallest magnitude of a quartic equation: 

8£ 2 7 ^[2 + (N- 2) 7co ] 2 - (1 - 7oo )(l + (N - l) 7co ) 3 = 0. 



The infinite- D radius and energy per atom are then given by 

eMi + (7v-i)7oor l/2 



l(7V-l)(7V-2)7^+2Ar 7o 



(95) 

(96) 
(97) 



N 2 (l-7oo)[l + (JV-l)7oo] ' 

where we use the superscript (QD) to distinguish similar quantities in the atomic and hard-sphere systems. 
Just as in the iV-electron atom analysis, we have used the fact that we are considering a minimum in 
the effective potential, which allows us to eliminate the coupling constant £ from the above expressions. 
Specifically, we have employed the conditions (20) and (21) at the minimum that 



din 



which leads to 



2 3/2 7oo[2+(A-2) 7co ] 



(98) 
(99) 



VT^-[1 + (N - l) 7 oo] 3 / 2 ' 
6.2 iV-electron Quantum Dot Normal Modes 

We can determine the elements of G by comparing the differential term in Eq. (30) with 7~( QD ) in Eq. (89) 
expanded to first order in 1/D as in Eq. (29). The non-zero elements of the G matrix are: 

G a = 1 

G„ = 2-£*£ = 2(l- 7oo )(l+7 00 )[l + (JV-l)7 00 ] 
(h = ~~T~~^~~rr~~ = 7oc(l - 7oc)[l + (-V - l) 7 oo], 



(100) 



7oo(l-7oo) 
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where the matrix elements have been evaluated at the infinite- D symmetric minimum and we have used 
the notation in Eq. (34). (See the discussion after Eq. (75) for an explanation of the factor of 2 in G g .) 
Likewise, the non-zero F matrix elements are: 

Fa = on _ X , 2 (8 + 2(5iV - 13) 7 co + (57V 2 - 217V + 24) 7 ^ - 3(N - 1)(N - 2)'&) 



F b 



"Too 

( 7oo -3)(2 + (7V-2) 7oo ) 



2(1 -7oo) 



2 



f - = 2(i- 7 ^ ( r + 7 ^-ih 0o) v 2 ( 2+ ^- 2 ^ 

Ff = (i- 7 oo) 2 (i+Tjv- ih^F 5 (101) 

F g = 2(1 - 7oo)3(i I {N - ihoo)2 (4 + 6(2N - 5) 7 ^ + (64 - 53N + UN^ + 
(-38 + 497V - 23iV 2 + 47V 3 ) 7oo 3 ) 

F " = (i - w ( itT« - ih~)» (3 + (5JV - 14h ~ + 111 - m + 2A,2h ~ > 

« - ^f-^^(2 +(W -2 h „), 
y 1 loo) 

where we have used Eqs. (D13) from Appendix D to evaluate the Gramian second-order derivatives. Again, 
since we are considering a minimum in the effective potential, we employed Eq. (99) to replace the coupling 
constant £ of Eq. (90) in favor of 7oo . Using the above equations for F and G along with Eqs. (42) we find: 

a = — 1 (8 + 2(5iV - 13) 7oo + (5iV 2 - 217V + 24) 7 ^ - 3(N - 1)(N - 2) 7 ^) 

^{l 7ooJ 

, _ 7co(7co-3)(2 + (Ar-2) 7co ) 



c = 



2(1 - 7 oc) 2 

AT _ 1 U/2 

(4 + 3(/V-2) 2 7co 2 ) 



-7oo(l+7oo)(l + (A r -l)7oo) 1/2 , ,, > 



2(1 -7oo) 



d = ^ 2(l + ^ ) 1 (l + ( ^ lh - )l/2 (2 + ( j V-6) 7co ) 
(1 - 7oo) 



e = 



(2 + 3(iV-2) 7oo ) (102) 



/ 
5 



2(1 -7co) 2 (l + (^V- 1)700)^ 
2^ 

(l-7co) 2 (l + (^-l)7oo) 1 / 2 
1 /. 



1 = 



(l-7oo) 2 (l + (^V- 


l)7oo) V 


(-30 + 20iV + 277V 2 


- 12iV 3 


Too 




2(l- 7oo ) 2 (l + (/V- 


- l)7oo) 


(-14-35iV + 457V 2 


- 14iV 3 


-4 7oo 2 (l +7eo 

t-l \2Z1 , T AT 


> i 



(4 + 6(2IV - 7) 7co + (16 - 17/V + 2iV 2 + 2iV 3 ) 7co 2 + 



(-1 - (TV - 6) 7co + (37V - 5) 7co 2 ) . 



Substituting these GF matrix elements in Eqs. (32), (55), (56), and (58), we now write down the quantum 
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dot normal-mode frequencies. The multiplicity- TV (iV — 3)/2 mode is a vector mode with frequency 

4 m) = Vg-zh + L. 



(103) 



The multiplicity- (N — 1) asymmetric stretch and bend frequencies take the form 



i- 



(104) 



and the two multiplicity- 1 symmetric stretch and bend frequencies take the form 



= \A?o ± Vm 2 - A, 



•o- 



(105) 



where r)\ and Ai are given in Eq. (57), and and A are given in Eq. (59). 

Having obtained the frequencies from Eqs. (103), (104), and (105) and the Lewis structure energy in Eq. 
(97), the energy through first order is given by Eq. (60). Although this is only a low order approximation, 
two-body studies suggest that when the coupling constant £ << 1 (i.e., the strongly interacting regime where 
the repulsive energy of the electrons dominates the confinement energy of the trap) low order dimensional 
perturbation may be very accurate[13]. Dimensional perturbation theory has a non-perturbative character 
in the sense that the leading-order term of DPT includes a contribution from the Coulomb potential. 

7 N hard-sphere particles in a trap: Atomic Bose-Einstein con- 
densates 

In this section we treat N hard spheres in a trap. An important application of this many-body system 
is N trapped ultra-cold Bose atoms (i.e., an inhomogeneous Bose-Einstein condensate) [17]. Dimensional 
perturbation techniques were introduced to inhomogeneous atomic BEC in Ref.[18] by scaling the non- 
linear mean-field (Gross-Pitacvskii) equation. The treatment in this section goes beyond the mean-field 
approximation as it includes correlation and uses a finite-range, as opposed to a zero-range, interatomic 
potential. 

The atoms of equal mass m are confined by a spherical harmonic trap with trapping frequency lu^o- 



We take the interaction potential to be a hard sphere of radius a (also the scattering length for the BEC 
system): 



I 0, Tij > a. 

In the two previous systems, many-electron systems, the dimensionally continued Laplacian is dimension 
dependent while the potential energy maintains the same form as it has at D = 3. In this system, we 




(106) 




(107) 
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dimcnsionally continue the interaction potential so that in the large- D limit the effective potential is diffcr- 
entiable, and thus allowing us to perform the dimensional perturbation analysis outlined in Sees. 3 and 4 
which requires taking first- and second-order derivatives of the effective potential. Specifically, we take the 
interaction to be 



V 



1-3/D 



1 — tanh 



1-3/D Vv^ 



I — a — —(a - a) 



D 



1 + (1-3/D) ^ 



*-3 2n 



n=l 



2™ 



(108) 



where D is the Cartesian dimensionality of space. This interaction becomes a hard-sphere of radius a in 
the physical, D = 3, limit. The other s constants (V , a, and {c„; Vn:0<ti<s - 3}) are parameters 
that allow us to fine-tune the large- D shape of the potential and optimize our results through Langmuir 
(first) order [16]. The simplest possibility could have as few as two parameters: V and c D , with a = a and 
the remaining c„ = 0; however, we can have any number of parameters for the most general and flexible 
potential. To see how the potential may be fine-tuned to optimize the results through Langmuir order, see 
Rcf. [16]. 

We use dimensionally scaled harmonic oscillator units similar to the TV-electron quantum dot system with 
the difference that we use D 2 in the dimensional scaling instead of of Eq. (86) . Ideally, one would use Q, in 
order to remove all of the dimension dependence from the centrifugal-like term in the kinetic energy, but in 
order to simplify the scaling of our dimensionally continued hard-sphere potential, we allow some dimension 
dependence in the centrifugal-like term (see [/( HS ) in (112) below). We regularize the large-D limit of the 
Schrodinger equation by using the following dimcnsionally scaled variables: 



E = £—E, H = ^-H, a 



V2D 2 



V — U V cv — 2 



, c = \/2D 2 a ho c , c„ = (V2D 2 a h0 ) 2n c n , 



where 



"Ho 



and iv h0 = D 3 uj ho 



(109) 



(110) 



are the dimensionally-scaled harmonic-oscillator length scale and dimensionally-scaled trap frequency, re- 
spectively. Substituting these scaled variables into the similarity transformed Schrodinger equation, Eq. 
(13), gives the following equation: 



(Y (HS) + c/ (HS) + y (HS) ) $ = 



(in) 



where 

T (HS) = A 5 2J2 



N 



d 2 

df 2 



EE 



Ijk - Hjlik 



d 2 



um = ^(s-i)[(2iy + i)a-i]r(' 



i=l 
N 



8f 2 



(112) 
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EE 



JV-l AT 



i=l i=i+l 



1-35 VV2 



3L_a-3*(a-a) 1 + (1 - 35) £ ^| 



«-3 g =2n ~ 



n=l 



where 5=1 /£) is the perturbation parameter and the interatomic separation is 



(113) 



As D becomes infinitely large, and 5 — > 0, the entire differential part of the kinetic energy as well as a portion 
of the interatomic and centrifugal-like potentials will drop out of the Hamiltonian. In the infinite-dimension 
limit, the particles behave as though they become infinitely heavy and rest at the bottom of the infinite-D 
effective potential, a potential which includes the confining potential and contributions from the centrifugal- 
like and interparticle potentials. The infinite-D energy becomes the minimum value of the effective potential 
(22). 



7.1 Hard-sphere infinite-D analysis 

The infinite- D (5 — > 0) effective potential in dimensionally scaled harmonic oscillator units is 

N / , _c-i , x JV-l N / 



i=l 



i=l j=i+l \ 



8-3 - =2n ' 



(114) 



8ff r 2 

As one can see from the double-sum term in Vg ff , the interaction potential becomes a soft sphere of radius 
approximately a and height 2V . The slope of the soft wall is determined by c OJ while, as discussed earlier, 
the remaining s — 3 parameters act to further refine the shape of the interaction potential[16]. 

Again choosing the totally symmetric configuration for which all f t and jij are equal to some and 
7oo, and using Eqs. (D12) from Appendix D in Eqs. (20) and (21), we find that the large- D radii and energy 
per atom are: 



E 



(HS) 



= [2(1 + (N - l) 7oo )]- 
1 + (N - 2) 7o 



-1/2 



1 



N (l-7co)(l + (A-l) 7co )8F^ + 2 r °° + 2 

where for simplicity of presentation we have defined the following: 



N — 1 



V [1 - tanh(9)] , 



s-3 



6 = C (f^y/T^Z- a) 1 + E C»^(l - loo)' 



(115) 
(116) 

(117) 



71=1 



The large- D direction cosine, 700, of the hyperangle between the infinite-dimensional radii is given by 
the negative solution of smallest magnitude of 

2(2+(7V-2) 7co ) 2 



V o c o Tsech 2 + 7o 



= 0, 



where 



(l-7co) 3 (l + (W-l)7co) 

1 + £ ((2n + l)c„e(l - 7co)" - 2nac„f 2 "- 1 (l - 7 co)"- 1/2 ) 



s-3 



n=l 



(118) 



(119) 
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7.2 Hard-sphere Normal Modes 

We can determine the elements of G by comparing the differential term in Eq. (30) with in Eq. (112) 

expanded to first order in 1/D as in Eq. (29). The non-zero elements of the G matrix are found to be: 

G a = 1 

G 9 = 2±-^ = 4(l- 7oo )(l+ 7oo )[l + (Ar-l) 7oo ] (120) 

G h = 7 °° ( r ^ 7 °° } = 2 7co (l- 7co )[l + (7V-l) 7co ], 
where the matrix elements have been evaluated at the infinite-!) symmetric minimum and we have used 
the notation in Eq. (34). (See the discussion after Eq. (75) for an explanation of the factor of 2 in G g .) 
Likewise, the non-zero F matrix elements are: 



F a = 1 + 



F b 



F, 



F h 
F, 



1 + (N - 2) 7o 



4r* (l-7oo)(l + (JV-l)7oo) 2 

s-3 

E 



] C °(7V-l)sech 2 e 



c (l - 7 oo)T 2 tanh6 + 



1 



CO 



2rooVl -7 
2naf^' 2 {\ - 7co )"- 1 [(^ - l)7co - n)\ ]c 



n I 

-2n-2/i „, \ra-lr 



*±±±?£-\l - 7 oo)"- 1/2 [(2n - l) 7oo - (2n + 1)] + 



yo5o -sech 2 e 



2 

s-3 



c (l - 7 oo)T 2 tanh6 + 



l + 7o 



2rooVl -7o 



+ 



2n+ 1 



E ^^(1 - 7co)- 1/2 [(2n + l) 7oo - (2n - 1)] + 



n=l 



-2naf 2 "~ 2 (l - 7oo )"- i [«7oo - (n - 1))] c„ 



1 + (AT - 2) 7o 



2f^(l- 7co ) 2 (l + (iV-l) 7oo ) 2 2 



H — scch 6 



-c rooT 2 tanh6 + 



+ — =L= + J2l [ (2« + 1) 2 ^(1 - 7co)- 1/2 - (2n) 2 aF 2 - 1 (l - 7C0)- 1 



2 V^^ ' ,tl 2 

Too 2 



2*1(1 " 7 oo) 2 (l + (7V-l)7oo) 2 
1 

2*1(1 ~ 7oo) 3 (l + (N — l) 7o o) 3 



(121) 



(l + 3(7V - 2) 7co + (13 - UN + 3iV 2 ) 7cc 2 + 



(N - 2) (4 - 37V + 7V 2 ) 7co 3 ) + ^sech 2 9 
- E f ^"f'V Hl - 7oo)- 3 / 2 + 



i^: T tanh9+ 2(i- 7co) 3/ 2 



+ 



n=l 



F2„(„ hof^(l- 7 oo) n - 2 

[3 + (57V - 14) 7oo + (11 - 9iV + 2N 2 ) loo 2 ] , 



4*^,(1 -7oo) 3 (l + (JV- l)7oo) J 
7^(27(^-2)7^) 

f^(l-7cx>) 3 (l + (iV-l)7oo) 3 



23 



where we have used Eqs. (D13) from Appendix D to evaluate the Gramian second derivatives. For trans- 
parency we leave the r M terms in Eq. (121) instead of using the explicit dependent form of f^. The 
above expressions for F and G are used to write the elements of GF [i.e., the scalar quantities a through 1 
in Eq. (42)]. We may then write the five normal- mode frequencies by means of Eqs. (55), (56), (58), and 
(32). For the multiplicity- AT(./V — 3)/2 mode, the frequency is given by 

w 2 HS) = \/g -2h + L. (122) 

For the two multiplicity- (N — 1) modes, the frequencies are of the form 

(HS) 



... - \A?i ± Vm 2 - Ai 

and the two multiplicity- 1 frequencies are of the form 



-.(HS) 

V 



\fm±V' 



A , 



(123) 



(124) 



where rji and Ai are given in Eq. (57), and rjo and Ao in Eq. (59). 

Because of the factors of 5 in the centrifugal-like and hard-sphere potentials and y( HS ) of Eq. 

(112)), there is a constant shift, v a , in the first-order energies: 



E 



<HS) 



E 



(HS) 



OO ^ 

E E K + ^,n„4 HS) + V ° 



2 } 



(125) 



where the n M are the vibrational quantum numbers of the normal modes of the same frequency u> 



(HS) 



and 



is the occupancy of the manifold of normal modes with vibrational quantum number and normal 

(HSl 

mode frequency uj^ (see the discussion after Eq. (60)). The FG matrix method gives the normal- mode 
frequencies of the infinite- 1? minimum structure. However, due to the explicit dimension-dependence (i.e., 
5 = 1/D) in [/( HS - ) and V^ us \ one must also take into account the explicit constant factors of order S that 
appear in the power expansion of Eq. (112) which are not the result of second derivatives of the effective 
potential. This contribution from the power series is the order-5 shift, v , in the energy in Eq. (125), and is 
given by: 



Vo = 



JV(JV + l)(l + (iV-2) 7co ) +3 y g N(N-l) 



4f 2 <J (l + (7V-l) 7oo )(l- 7(X) ) 



's-3 



(a -a) £c n e(l- 



V7l=l 




(126) 



where & is given by Eq. (117). 



8 Summary 

DPT has advantages over traditional methods that cither neglect the interaction at lowest order (con- 
ventional perturbation theory) or the kinetic energy at lowest order (e.g., the large- N Thomas- Fermi 
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approximation[14, 18]). The infinite-dimension limit (leading order) of DPT results in an effective po- 
tential that keeps contributions from the confinement, interaction, and kinetic terms of the Hamiltonian. 
The kinetic contribution that remains in the large- D limit is a centrifugal- like (1/r 2 ) term that allows the 
leading order of DPT to satisfy the minimum uncertainty principle. Correlation is included at all orders, 
including the large-dimension leading order of the 1/D expansion. In the infinite-dimension limit, we choose 
a completely symmetric configuration in which all of the particles are localized and equivalent. The first- 
order energy correction corresponds to small oscillations about this high-Z? symmetric structure, with the 
vibrational frequencies determined by the Wilson FG matrix method. There are five distinct roots of the 
GF matrix, belonging to three different irreducible representations^, 9] of the symmetric group of order TV, 
Sn- These five distinct roots are directly related to the normal mode frequencies of the many-body system. 
The two multiplicity- (TV — 1) roots, which are a mixture of asymmetric stretching and bending motions, 
are designated by 1" and 1 + , respectively. The two multiplicity- 1 roots, which are a mixture of symmetric 
stretching and bending motions, are designated by 0~ and 0+, respectively. The multiplicity- TV(TV — 3)/2 
root, designated by 2, is a purely angular mode. 

Analytical results through first order are obtained for a spherical system of identical particles with a 
general confining potential and general two-body interaction potential. Many-body DPT is applied to two TV- 
electron systems: the TV-electron atom originally studied in Loeser's seminal paper [1] where the confinement 
is supplied by the attraction of the nucleus, and the TV-electron quantum dot where the confinement is 
supplied by an external harmonic potential. We also consider the inhomogeneous bose condensate whose 
TV atoms interact via a hard-core potential. As in the TV-electron quantum dot, the condensate atoms are 
confined by an external harmonic trap. Unlike either TV-electron system, however, the hard sphere potential 
has explicit dimension dependence. Both the atom and quantum dot interact via a repulsive Coulomb 
potential that maintains its three-dimensional form as D varies (an alternative generalization of the three- 
dimensional Coulomb potential might be l/r D ~ 2 ). In order to make the hard sphere amenable to DPT 
analysis, we dimensionally continue the hard-sphere potential so that it is differentiablc in the infinite- D 
limit and becomes a hard-sphere with radius equal to the scattering length in the physical D = 3 limit. This 
dimensional continuation results in a soft-sphere interatomic potential at large D. 

While the treatment of these physical systems has been in many respects quite general, there remain 
many avenues for further development of this many-body DPT formalism. One extension might be to relax 
the spherical symmetry constraint by allowing confining potentials with axial symmetry. This is a particu- 
larly important generalization for describing many experiments, especially Bose-Einstein condensates, which 
in practice are characterized by axial symmetry. With axial symmetry, one has TV additional coordinates due 
to the z-components of the TV particles, but the many-body formalism presented in this paper readily accom- 
modates this generalization. Instead of generalized spherical coordinates with a D-dimensional radius, one 
uses generalized cylindrical coordinates, which consists of a (D — l)-dimensional sphere plus a ^-coordinate. 
Another extension involves the choice of confining potential. This paper focuses mainly on harmonic traps 
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because this is the most prevalent form of confinement in experiments today, but there is nothing to prevent 
one from also applying these methods to more general external confining potentials. Yet another extension 
involves the calculation of higher-order perturbation terms. The leading-order DPT wavefunction is a prod- 
uct of one-dimensional harmonic oscillator wavefunctions with frequencies given in the first-order term of 
the energy series (from the FG method). From this wavefunction, one can then use ordinary perturbation 
theory to extend these results to include corrections beyond first order. This is also an important extension 
as it should improve DPT's numerical results for many-particle systems. Even though the 1/D perturbation 
series is asymptotic, with appropriate treatment higher order terms can be used to improve the accuracy of 
the low order approximation. 

Numerical results are given in Refs. [1] and [16] for the ./V-electron atom and TV-atom condensate, 
respectively. Numerical analysis for a two-electron quantum dot in Ref. [13] suggests that the low-order 
analytical many-body results given in Sec. 6 will be quite accurate for typical many-electron quantum dots 
whose length scales are in the mesoscopic regime (i.e., the strongly interacting or strongly confined regime). 
Moreover, because N is a parameter in our formalism, the challenge of calculating the physical properties of 
a given system with DPT does not increase as one adds more particles. The methods detailed in Sees 2-4 
can be applied to any variety of many-body systems, and we hope that the systems described in Sees 5-7 
are sufficiently varied to encourage further interest in this novel many-body approach. 
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A Wilson FG matrix method 

In this appendix, we derive the Wilson FG matrix method [2] which is at the heart of our obtaining the normal 
mode frequencies and, thereby, the first-order energy correction. The derivation involves a transformation to 
the set of coordinates called normal-coordinates in which both the differential term and the potential term 
of Eq. (30) are diagonal. We begin by defining a symmetric transformation, A, that transforms from the 
vector y', defined by Eq. (26), to z! (z' — Ay'). A is an active transformation and it has the property that 
it diagonalizes the symmetric G matrix to unity. That is, A satisfies 



A T GA = I, 



(Al) 



which can also be written as 



G = A-^A^f, 



(A2) 
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where we have used the property that A is symmetric (A = A T ) in the derivation of Eq. (A2). Hi then 
becomes 

H 1 -\d?,d- z , + VtA-TFA-V. (A3) 

Next we focus our attention on the potential term, the term involving the matrix F. We introduce 
another transformation, U (q' = Uz'), that diagonalizes the potential term while simultaneously leaving 
the differential term unchanged. This orthogonal transformation (U T U = I, where I is the identity matrix) 
leaves the differential term in the same form as in Eq. (A3), and the potential term becomes 

U(A- 1 ) T FA- 1 U T = A, (A4) 

where A is a diagonal matrix. That is, 

^i--^^ + iq' T Aq'. (A5) 

The eigenvalue equation corresponding to Eq. (A4) is 

(U(A- 1 ) T FA- 1 U T ) q' = Aq', (A6) 

where q' is the eigenvector with eigenvalue A. Equation (A4) is a matrix eigenvalue equation expressed in 
the basis of the normal coordinates. We can change the basis back to the original internal displacement 
coordinates with the passive similarity transformation q' = A _1 U T q', under which Eq. (A6) now reads 

U(A- 1 ) T Fg' = XXJAq'. (A7) 

Multiplying on the left by U T followed by A -1 gives 

A- 1 (A- 1 ) T Fy' = Ag'. (A8) 

Inserting Eq. (A2), gives the eigenvalue equation for the normal mode coordinate q': 

GFq' = \q'. (A9) 

Thus, the FG matrix method consists of finding the roots of the characteristic polynomial in A: 

det(AI — GF) = 0, (A10) 

which is carried out for a general A-body system in Sec. 4 and applied to the systems of Sees. 5-7. 



B Spectral graph theory and the symmetry of the Q matrices 

In this appendix, we introduce the relevant aspects of graph theory in order to derive several quantities used 
throughout the main text and to show the connection of graphs to the large- 1? symmetric configuration of 
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N particles. A more complete treatment of the mathematics of graph theory with an emphasis on chemical 
physics applications may be found in Ref. [10]. 

A graph Q = (V,E) consists of a finite set of vertices {v\,V2, ■ ■ ■ ,v n }eV(G) and a set of not necessarily 
distinct unordered pairs {(vi, Vj), V(l < i < j < n)}eE(G) connecting the vertices to form the corresponding 
edges. A graph is complete on n vertices, denoted by K n (Q), if every distinct pair of vertices is connected by 
an edge; K n is sometimes called a simplex of n points. Notice that K n contains n(n — l)/2 edges. Our large- 
dimension symmetric minimum configuration of N atoms forms a simplex on TV vertices in D-dimcnsional 
space. For example, the N — 4 simplex is a hypertetrahcdron. The line graph L(Q) of a graph Q is the graph 
whose vertices correspond to the edges of Q with two vertices being adjacent if and only if the corresponding 
edges in Q have a vertex in common. 

Derivation of Eq. (42): The Q matrices display a high degree of symmetry in the large-dimension 
configuration, but it is not trivial to write down the elements of GF in terms of its constituents G and F 
as we have done in Eq. (42). We now derive this result in a two-step process. The first step is to write the 
(N + M)-square matrices G and F, where M = N(N — l)/2, as block matrices of the following form: 

Q=( QlQi V (Bl) 
V Q 3 Q 4 J 

where the blocks, each of which will be discussed in turn, have the same indical structure as the blocks of 
yy T in Eq. (33). 

From Eq. (34), the (N x N) Qi block is composed of two distinct elements: the diagonal elements 
Q a = Fa and the off-diagonal elements Qb = Qij (i =/= j)- Qi can be written in terms of the more basic 
matrices Ijv (the (N x N) identity matrix) and Jn (the (N x N) matrix consisting of all ones): 

Qi = (Q a - Q b )l N + QbJ N . (B2) 

The matrix Jjv contains all ones, so for the diagonal parts of Qi the Qb's cancel and one is correctly left 
with Q a . The off-diagonal parts of Qi correctly become Qb since the off-diagonal of Ijv is zero. Specifically, 
taking the Q matrix to be F, we find 

Fi = (F a - F b )I N + F b 3 N . (B3) 

One can write down the Gi block in similar fashion noting that the off-diagonal entries of Gi are zero: 

Gi = G a l N . (B4) 

The (N x M) Q2 block according to Eq. (34) is also composed of two distinct elements: the "incident" 
elements Q e = Qi^j (i 7^ j) and the "non-incident" elements Qf — Qi.jk (i j =/= k). The elements F e 
are termed "incident" because the repeated index i means that the vertex i is incident with edge (ij), while 
for Qf the vertex designated i is not incident with the edge (jk) since i =/= j =/= k. To write Q2 in terms of 
basic matrices, we need the (N x M) vertex-edge matrix R, defined as Rijk — &ij + <5ife, which equals one 



28 



when vertex i is incident with edge (jk) (i.e., when i = j or i = k) and zero otherwise. Thus, the following 
equation accurately describes the Q2 block: 

Q2 = {Qe-Qf)n + Q f 3 NM . (B5) 

Specifically, letting the Q matrix be F, we find 

F 2 = (F e -F f )TL + F f J NM . (B6) 

And since there is no mixing of radial and angular derivatives in the Schrodinger equation (13), we have 

G 2 = G 3 = 0. (B7) 

By definition, the Hessian matrix F is symmetric; thus, from Eq. (B6) we find 

F 3 = (F e - F f )R T + F f 3 MN , (B8) 

where Jmn = Jjvm- However, for a more general Q-matrix (e.g., GF which, unlike F, is not symmetric), 
an analysis analogous to that leading to Eq. (B5) yields 

Q 3 = (Q c - Qd)R T + Q d JMN, (B9) 

where Q c = Qij t i (i j) and Qd = Qjk,i (i =/= j =/= k). The Q3 block can be interpreted in the same was as 
Q2 where the "incident" elements are now Q c and the "non-incident" elements are Qd 

For the fourth block of Q, we need to introduce the adjacency matrix of the line graph of the simplex, i.e., 
the quantity R T R — 21m- A more intuitive way to interpret this quantity is to think of it as an edge-edge 
matrix that is unity when two edges are adjacent and zero otherwise. For the purposes of constructing the Q4 
matrix, which has indices of the form (ij, kl), it helps to then think of the index (ij) as an edge of the simplex 
connecting vertices i and j, and likewise for edge (kl). Two edges are adjacent if they share a vertex. Q4 is 
comprised of three distinct elements: the diagonal elements Q g = Qij,ij (i =/= j), the adjacent-edge elements 
Qh = Qij.jk (i =/= j k), and the non-adjacent-edge elements Q L = Qij^ki (i ^ j ^ ^ /), which have no 
repeated indices and, therefore, no adjacent edges. The matrix R T R, whose (ij, kl) th element is given by 
hk(JN)ji + Iu(JN)jk + Ijk(JN)u + Iji(JN)ik, gives the adjacent edges, but it double counts the diagonal; 
thus, the term Q/j(R T R — 21m) of Q4 accounts for the adjacent edges. To get the diagonal elements of Q4, 
we need the term Q 9 1m- And for the non-adjacent edges we need Q l [Jm — (R T R — 1m)]- Putting these 
terms together gives: 

Q 4 = (Q g - 2Q h + QJI M + (Qh - QJR T R + QJm- (BIO) 

Specifically, one has 

F 4 - (F g - 2F h + FJI M + (F h - F t )R T R + F t J M , (Bll) 
and likewise for G4, noting that G L = 0: 

G 4 = (G g - 2G h )l M + G h R T K. (B12) 
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The second step to derive Eq. (42) is to find the blocks of GF by matrix multiplication of G and F, 
whose blocks are given above. That is, we multiply the following 

/ Gi \ / Fi F 2 \ 

G = F = . (B13) 

\ G 4 J \ F 3 F 4 J 

Then by analyzing each block of GF we can write down its elements in terms of the elements of G and F. 
The first block of GF gives a and b of Eq. (42) and is given by 

(GF)! = GiFi = (G a F a - G a F h )l N + G a F b 3 N . (B14) 

Taking the diagonal part of (GF)i gives 

a = (GF) a = {GF)i,i = G a F a , (B15) 

while the off-diagonal (i ^ j) gives 

b = (GF)b = (GF)ij = G a F b . (B16) 

The second block gives e and / of Eq. (42) and is given by 

(GF) 2 = G!F 2 = (G a F e - G a F f )K + G a F f 3 NM . (B17) 

Noting that Ri^j is unity when i ^ j and that all elements of Jnm are unity, we find the following from 
(GF) 2 when i ^ j: 

e = (GF) e = (GF) iti j — G a F e . (B18) 

For i ^ j ^ k, Rijk = 0; and so, 

/ = (GF) } = (GF)ijk = G a Ff. (B19) 
When deriving (GF) 3 and (GF) 4 below, the following relations prove useful: 

RR T = Jn + (N — 2)Ijv R t J w = 2J M jv RJmjv = (N - 1)3 N (B20) 
RJm = {N — 1)3 nm R t Jatm = 2Jm • (B21) 

Equations c and d of Eq. (42) can be found from the third block of GF, which is given by 

(GF) 3 = G 4 F 3 = [G g F e - G g Ff + (N - 4)(G h F e - G h F f )} R T 

+ [G g F f + 2G h F e + 2(N - 3)G h F f ] 3 MN . (B22) 

For the (ij, i) elements (i ^ j), R T is unity; thus, to get (GF) C one adds both bracketed terms in Eq. (B22): 

c = (GF) C = {GF) ij}i = G g F e + (N — 2)G h F e + (N - 2)G h F f . (B23) 

For the (jk, i) elements (i ^ j ^ k), R T is zero; thus, {GF)d is simply given by the second bracketed term 
in Eq. (B22): 

d = {GF) d = {GF) jkti = GgF f + 2G h F e + 2(N - 3)G h F f . (B24) 
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Equations g, h, and i of Eq. (42) can be found from the fourth block of GF, which is given by 

(GF) 4 = G 4 F 4 = \{G a - 2G h )(F g - 2F h + F)] 1 M + [G g F L + AG h F h + 2(N - A)G h F L ] 3 M 

+ [(N - 6)G h F h — (N — h)G h F b + G g (F h - F L ) + G h F g ] R T R. (B25) 

As noted earlier, R T R double counts the diagonal elements (i.e., elements where i 7^ j). Thus, the 

(GF)ij^j elements are given by the sum of the first two bracketed terms in Eq. (B25) plus twice the third 
term: 



g = {GF) lhlJ = (GF)ij t ij = G g F g + 2(N - 2)G h F h . 



(B26) 



For the (ij,jk) elements (i ^ j ^ k), we note that the R T R elements are unity and the elements of Ijv are 
zero. Thus, to get (GF)ijjk we take the sum of the second and third bracketed term in Eq. (B25): 



h = (GF) h = (GF)ijjk - G g F h + G h F g + (N — 2)G h F h + (N - 3)G h F t . 



(B27) 



Finally, the (ij, kl) elements (i ^ j ^ k ^ I) of R T R and Ijv are zero, and thus the last equation in Eq. (42) 
is simply given by the second bracketed term in Eq. (B25): 



1 = (GF), = (GF) ijM = G g F L + AG h F h + 2(N - A)G h F L . 



(B28) 



Derivation of Eq. (49): A graph with a given structure corresponds to a spectrum. That is, a graph 
may be represented by a matrix whereby its spectrum, or eigenvalues, may be calculated. Pg(X) denotes the 
characteristic polynomial of the adjacency matrix A of graph Q: 



Pg(\) = det(AI — A). 



(B29) 



The element of the adjacency matrix is the number of edges connecting vertices i and j of the graph Q. 
The spectrum of Q is found by solving Pg(X) = 0. 

To prove the relationships involving determinants in this appendix and the next, we quote the following 
result from p. 72 of Ref. [10] for the complete graph K n : 



PkJX) = (A-n + l)(A + l) 



71-1 



(B30) 



The adjacency matrix of a simplex of n points, K n , is an n x n matrix of ones except along the diagonal 
which contains all zeros. That is, Pk„ (A) takes the following form: 

/ A -1 ••• -1 \ 

-1 A '•• : 

: *•. '•• -1 

V -1 ••• -1 A / 



P Kn (A) = dct 



(B31) 
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Using Eq. (B30) one can easily derive Eq. (49) of Sec. 4 , which states 



vl 



M 



det(XYEZ) = det y " 
where the integer M is defined in Eq. (35). First we note that 
det 



= t A, - 1 (t + Nu)v 



M-N 



(B32) 



Vl 



det (vIm) det I -I N + —Jn 



A I 



V V 



(B33) 



Then, det(^I]v + ^jJjv) can be written as an (N x N) matrix of the same form as Pjf n (A) in Eq. (B31): 



f t+u u _ u ^ 



, . t u . 

det(— Ijv H — Jat) = det 

v v 



u t+u 

V V 



(-- 

V V 



N 



det 



( -l-i- -l 

-1 -l-i 
V -l 



-1 \ 



-1 -1 - i- J 

u / 



(B34) 

so that, noting det(wlAf) = v M , and setting A = — 1 — £ and n = TV in Eq. (B30) leads directly to 
det(XYEZ) = t N -\t + Nu)v M - N . 



C Gaussian Reduction of Equation (44) 

The characteristic equation for the GF eigenvalues is det(E) = 0, where E is given by Eq. (44). The 
Gaussian elimination to reduce E of Eq. (44) to the form XYEZ of Eq. (45) allows us to find analytical 
expressions for the normal-mode frequencies. The goal is to transform E to a lower-triangular matrix, whose 
determinant we can compute exactly with Eq. (49), while leaving the characteristic determinant invariant 
by imposing that the transforming matrices have unit determinant. The elimination process consists of three 
steps. 

Step 1: We first define Z: 

(Ijv wJnm \ 
(CI) 
I M J 

whose determinant is unity, which we now show. It is known from matrix theory (for example, see Ref. [11]) 
that the determinant of an (N + M)-square matrix with an (N x M) zero matrix (Ojv,m) m the upper right 
block or an (M x N) zero matrix (Om,n) in the lower left block is the product of the determinant of the 
diagonal block of matrices. That is, 





OjVM 




f A 




; 


= det 




I 








[ Omn 





det =det = det(A) det(B), (C2) 

\C B J \ Omn B J 

where A and B are square matrices. For Z, we have det(A) = det(Ijv) = 1 and det(B) = det(lAf) = 1; thus, 
det(Z) = 1. 
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Multiplying E on the right by Z and using the multiplication rules above leads to 

E (A - a + b)I N - bJ N (f-e)-R+[w(\-a-(N-l)b)-f}J NM 

(d - c)R T - dJ NM (\-g + 2h-L)lM-(h-L)-R T ~R+[(2(d-c)-Nd)w-i}J M 

We choose w equal to //(A — a — (N — l)b), such that the upper-right Jnm term vanishes. This results in: 

(A -a + b)I N -bJ N (/-e)R 
EZ = ( V U ; | , (C4) 

(d - c)R T - dJ NM v! M - (h - i)R T R + C,3 M 



where 



A- a - (N- l)b y 

v = A - g + 2h - l. (C6) 

Step 2: Next multiply on the left by Y, which by Eq. (C2) has unit determinant: 

V-f °V (C7) 

y xR T + yJ N M T Im / 

where we impose 

2(d-c) - Nd 



_ fr-z _ (_ 1_ 



A -a- (W- 1)6 



to eliminate the lower-right R T R and J m terms from EM. This results in 



where 



(08) 



YEZ = | + C/-.)H , (C9) 

£R T + mJ NM i>I M 



Z = ^ — -(A- a + 6) + (d- c) 
/ - e 

to = y(A — a — (N — 1)6) - 26-^— - (CIO) 

/ - e 

-[/(2(d - c) - Mf) - i(A - a - (JV - 1)6)] - 2b^—^- - d. 



Step 3: The final step is to eliminate the upper right term from YEZ by multiplying YEZ by X, which 
by Eq. (C2) has unit determinant: 

I I N zR \ 

X^ , (Cll) 

\ I M J 

where z = (e — f)/v. XYEZ results in the desired form of Eq. (45): 

, ±i n + -Jn o , 

XYEZ = ( v v | , (C12) 

ZR T + mJ NM v! M 
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where 



* = {X-a + b)v + {N-2)(t-h)(X-a + b) + {N-2)(d-c){e-f) (C13) 
u = -kv- (h- t)(A- a + b) + {d- c)(e- /) 

+ [/(2(d - c) - JVd) - t(A - a - {N - 1)6) + 4b(h - t) - 2d(e - /)] . (C14) 

As the determinants of Z, Y, and X are unity, the determinant of XYEZ and, hence, the characteristic 
determinant det(E) may be calculated using Eq. (49), which was derived in the previous appendix. 



D Gramian Determinants 



The Gramian determinant [19] is defined as: 



det(7y), 



(Dl) 



where 7 y = • Vj/riVj, the angle cosines between the particle radii r 4 and rj, represents the elements of an 
N x N matrix. A related quantity used in the main text is the principle minor of the Gramian, , defined 
as the determinant of the 7^ matrix with the a th row and column removed. 

A most challenging part of calculating the large- .D minimum and the F matrix elements in the systems 
discussed in this paper is handling the Gramian determinants and their derivatives. What makes these 
calculations feasible is the very high symmetry of the infinite-dimension, symmetric minimum. We make use 
of Eq. (B30) to obtain the Gramian determinant and its derivatives at the infinite-dimension, symmetric 
minimum. We will now demonstrate how this is done and summarize the results. 

In general 7^ = 1 and at the infinite- 1? symmetric minimum all of the remaining direction cosines are 
equal, 7^ = 700. Hence, the Gramian determinant is an (N x N) matrix of the form 



det 



/ 1 

V 7oc 



7o 



7o 



7o 



(D2) 



1 / 



which can be written in the same form as P^ re (A) (Eq. (B31) of Appendix B): 

( -l/7oo -1 ••• -1 \ 
-l/7oo '•• : 



(-7oo) W det 



-1 



V -1 



-1 



-1 



(D3) 



-l/7oo / 



Setting A = — 1/ 700 and n = N in Eq. (B30) gives 



[l + (JV-l) 7oo ](l- 7oo ) 



JV-l 



(D4) 
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The principle minor evaluated at the infinite- D symmetric minimum 1^"' is simply related to the corre- 

oo 

sponding Gramian determinant (Eq. (D4)) by N — > N — 1. 

To calculate the infinite- D symmetric minimum of the Gramian derivatives, we expand T in terms of its 
cofactors. The cofactor, denoted by dj, of the element 7^ in T is (— multiplied by the determinant of 
the matrix obtained by deleting the i th row and j th column of T. We may then write T as 



N 

r = 7, j Cij . 

Then the partial derivative of T in terms of the cofactor is 

dT -2C 

From this equation, the partial derivative of T evaluated at the infinite- D symmetric minimum is 

dT 



(D5) 



(D6) 



(D7) 



where we have defined the following determinant of an (N — 1) x (N — 1) matrix: 

^ -~v -~v -y . . . >y ^ 



C^- 1 ) = det 



Too 


Too 


Too 


Too 


"Too 


1 


Too 




"Too 


Too 


1 





(D8) 



y 7oo Too 1 

and the superscript (N —1) simply indicates the size of the matrix. From this matrix one can show that the 
following recursion relation holds: 



(iV-2) 



e^TJr ~ -(N-2)Cg 



-2) 



or equivalcntly 



- 7co (r 

(JV) 



(JV-1) 



-(iV-l)C^- 1 ) , 



(D9) 
(D10) 

. From the recursion 



where the (N) superscript in the notation T again refers to the size of the matrix T 

00 . 

relation (Eq. (D10)), one can easily prove by induction the conjecture that = 7oo(l — 7oo) JV ~ 2 and so 

dT 



^Too ( 1 Too ) 



7V-2 



(Dll) 



The derivative of the principle minor evaluated at the infinite- D symmetric minimum is simply related to 
the corresponding Gramian determinant derivative (Eq. (Dll)) by N — > N — 1. 

To summarize the above results, the following expressions are needed when calculating the minimum of 
the effective potential (Eqs. (20) and (21)): 



= -2 7oo (l- 700)^2 

[l + (^-l)7oo](l-7co) Ar - 1 



dr 



(a) 



= -2 7 cx>(l-7cx>) JV - 3 

= [l + (7V-2) 7o o](l-7co) Ar - 2 . 



(D12) 
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And when evaluating the F matrix elements at the infinite- Z) symmetric minimum, the following six 
second-order derivatives of the Gramian determinants are needed: 



3-T 



9 2 r<°» 



o 



o 92r(Q) 



\N-3 



oo 



-2{l- loo ) N -i{l + {N-Z) loo ) =2 7oo (l-7cx>) ; 

nj ijk qq 

(D13) 

:-2(l- 7co )^- 4 (l + (7V-4) 7co ) ggg^ = 2 7oo (l- 7oo ) JV " 4 . 
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